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ABSTRACT 


Studies  of  the  atmospheric  boundary  layer  always  require  values  for  the  dynamic  and  thermodynamic 
properties  of  the  fluids  present:  namely,  air,  water  vapor,  and  other  trace  gases.  Studies  of  the  marine 
boundary  layer  frequently  require  similar  properties  for  the  ocean,  in  addition.  This  handbook  collects 
functions  for  calculating  the  values  for  these  dynamic  and  thermodynamic  variables.  For  example,  it 
includes  equations  for  air  density,  water  vapor  density,  and  seawater  density.  It  also  describes  various 
humidity  variables  and  shows  how  to  convert  among  these.  It  includes  equations  for  the  specific  heats  of 
dry  air,  water  vapor,  moist  air,  seawater,  and  ice  and  for  the  latent  heats  of  water  in  its  various  phases. 
The  handbook  goes  on  to  show  how  to  calculate  the  molecular  transport  properties  such  as  the  thermal 
conductivity  of  air  and  the  diffusivities  of  heat,  water  vapor,  and  several  environmentally  important  gases 
in  air.  This  discussion  naturally  includes  values  for  the  molecular  Prandtl  and  Schmidt  numbers.  Finally, 
the  handbook  concludes  with  summaries  of  several  wind  force  scales  that  descriptively  categorize  winds 
on  land  and  sea:  the  Beaufort  Scale,  the  Saffir- Simpson  Scale  for  hurricanes,  and  the  Fujita  Scale  for 
tornadoes. 
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1  INTRODUCTION 

Studies  of  geophysical  boundary  layers  always  require  kinematic  and 
thermodynamic  constants  for  the  fluids  involved.  The  most  obvious  examples  are 
for  the  fluid  densities:  air,  pure  water,  seawater,  and  water  vapor.  Boundary-layer 
studies  frequently  involve  exchanging  properties  across  interfaces;  consequently, 
molecular  properties  like  the  kinematic  viscosity,  thermal  conductivity,  water 
vapor  diffusivity,  and  surface  tension  are  also  necessary.  In  turn,  the  ratios  of  the 
molecular  transport  quantities — the  Prandtl  and  Schmidt  numbers — are  recurring 
variables. 

Here,  I  summarize  the  values  and  functions  for  these  and  other  quantities  that 
I  have  found  to  be  accurate  and  useful  for  my  work  in  boundary-layer  meteorol¬ 
ogy.  This  handbook  is  admittedly  not  all-encompassing.  For  example,  I  tend  to 
focus  on  the  quantities  and  the  range  of  conditions  for  studies  of  air-land,  air- 
sea,  and  air-sea-ice  interaction.  I  also  describe  some  microphysical  quantities 
used  in  studies  of  aqueous  solution  droplets,  like  cloud  droplets  and  sea  spray.  I 
tend  to  simply  describe  ways  to  calculate  the  quantities  of  interest  without  also 
explaining  why  you  might  want  this  quantity.  Therefore,  check  any  good  book  on 
atmospheric  thermodynamics,  such  as  Iribame  and  Godson  (1981),  Pruppacher 
and  Klett  (1997),  or  Bohren  and  Albrecht  (1998),  for  academic  discussions  of  the 
thermodynamic  quantities. 
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2  AIR  DENSITY 

Dry  air  obeys  the  ideal  gas  law, 


=  m^p 

RT 


(2.1) 


where  pd 
Ma 
P 
R 
T 


density  of  dry  air 
molecular  weight  of  air 
barometric  pressure 
universal  gas  constant 
air  temperature  (K). 


At  standard  temperature  (T  =  T0  =  273.15  K)  and  pressure 
(P  =  P0  =  1013.25  mb),  (2.1)  gives  pdo  =  1.2922  kg  nf3.  Therefore,  (2.1)  also 
gives 


Pd 


\ 

'  p  ^ 

f  Tn  ^ 

r  p  ' 

=  1.2922 

0 

/ 

It  J 

(2.2) 


where  pd  is  in  kg  m  3  when  P  is  in  millibars  and  T  is  in  kelvins. 
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WATER  VAPOR  DENSITY 


Water  vapor  in  the  atmosphere  also  obeys  the  ideal  gas  law, 


Pv  = 


Mwe 

RT 


(3-D 


where  pv  =  water  vapor  density 

Mw  =  molecular  weight  of  water 

e  =  partial  pressure  of  the  water  vapor. 

To  be  rigorous,  we  need  to  recognize  that,  in  (2.1),  the  partial  pressure  of  dry 
air  is  not  the  barometric  pressure  P  but  rather  is  approximately  P  -  e.  That  is,  in 
rigorous  usage,  (2.1)  should  be 


Pd 


Ma  (P-e) 


RT 


(3.2) 


By  rearranging  (3.1)  and  (3.2),  we  can  write  the  following  expression  for  the 
barometric  pressure: 


P  =  (P-e)  +  e 


RTPd  !  RTPv 
Ma  Mw 


or 


P  = 


RTPil 

Ma 


1  + 


VMw 


(3-3) 


(3.4) 


Here, 


Pa  =  Pd  +  Pv  (3-5) 

is  the  density  of  moist  air;  and  we  recognize  pv/pa  as  the  specific  humidity,  Q 
(more  on  this  soon). 

Equation  (3.4)  rearranges  to 
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Pa  =  - 7 - 1 - r  ,  (3.6) 

RT(l  +  0.608Q) 

where  Ma/Mw  -  1  =  0.608.  Equation  (3.6)  implies  that  (2.1)  is  inaccurate  by  the 
factor  (1  +  0.608Q)  if  we  want  the  total  air  density.  Since  for  normal  atmospheric 
conditions  Q  is  seldom  larger  than  0.035  kg  kg-1,  the  term  in  parentheses  in  (3.6) 
is  always  between  1.000  and  1.022.  Therefore,  (2.1)  may  be  accurate  enough  for 
many  purposes. 

Nevertheless,  we  often  rewrite  (3.6)  as  (e.g.,  Lumley  and  Panofsky  1964,  p. 
214) 


Pa  = 


MaP 

RTV 


(3.7) 


to  preserve  the  form  of  (2.1)  while  retaining  the  accuracy  of  (3.6)  by  defining  the 
virtual  temperature 


Tv  =  T(1  +  0.608Q)  . 


(3.8) 


In  effect,  Tv  is  the  temperature  that  dry  air  must  have  to  produce  the  same  density 
as  moist  air  at  the  given  barometric  pressure. 
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4  WATER  VAPOR  VARIABLES 

Many  types  of  instruments  are  available  for  measuring  atmospheric  water 
vapor.  Some  actually  measure  the  water  vapor  density  (the  absolute  humidity), 
but  many  measure  derivative  quantities  such  as  the  mixing  ratio,  the  dew-point 
temperature,  or  the  wet-bulb  temperature.  Hence,  analyses  and  data  reporting 
often  require  converting  among  the  different  water  vapor  variables. 
Schwerdtfeger  (1976)  and  Pruppacher  and  Klett  (1997,  p.  106f.),  among  many 
others,  give  good  summaries  of  water  vapor  variables. 


4.1  Vapor  Pressure 

If  there  is  a  fundamental  water  vapor  quantity,  I’d  say  it  is  the  vapor  pressure 
e.  And  for  computational  purposes,  the  saturation  vapor  pressure  esat  is  usually 
the  fundamental  variable. 

I  use  Buck’s  (1981)  three  equations  for  the  saturation  vapor  pressure  (cf. 
Brock  and  Richardson  2001,  p.  86ff).  For  vapor  in  saturation  with  a  planar 
surface  of  pure  water  at  temperature  T  between  -20°  and  50°C,  Buck  gives 


esat(T)  =  6.1121(1.0007  +  3.46xlO“6P)exp 


r  17.502T  A 


240.97  +  T 


(4.1) 


where  esat  is  in  millibars  when  P,  the  barometric  pressure,  is  also  in  millibars. 

For  saturation  over  water  at  much  lower  temperatures,  -40°  <  T  <  0°C,  Buck 
(1981)  gives 


esat(T)  =  6.1121(1.0007 


+ 


3.46xl0“6 


P) 


r 

exp 

v 


17.966T 

247.15+Ty 


(4.2) 


Use  this  relation,  for  example,  to  compute  the  saturation  vapor  pressure  in  clouds 
composed  of  deeply  supercooled  water  droplets. 

Finally,  if  the  vapor  is  in  equilibrium  with  a  surface  of  pure  ice  (or  snow,  for 
instance),  Buck  (1981)  recommends  the  following  equation  for  -50°  <  T  <  0°C: 


22.452T 


V 


J 


esat(T)  =  6.1115(1.0003  +  4.18xlO“6P)  exp 


272.55  +  T 


(4.3) 
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As  an  alternative  for  the  saturation  vapor  pressure  over  ice,  Murphy  and 
Koop  (2005)  give  a  more  complex  expression  that  extends  over  a  wider 
temperature  range,  -165. 15°  <  T  <  0°C; 


esat(T)  =  0.0l[l  +  10~5P (4.923  -  0.0325T  +  5.84xl(T5  T2)] 

5723.265 


>exp 


9.550426  - 


T 


+  3.530681n(T)  -  0.00728332T 


1  -(4.4) 


This  gives  esat  in  millibars  for  P  in  millibars  and  T  in  kelvins.  Equations  (4.3)  and 
(4.4)  differ  insignificantly  over  their  common  range.  I  therefore  prefer  (4.3) 
because  it  is  mathematically  simpler.  Use  (4.4),  however,  for  temperatures  below 
-50°C. 

Figure  1  shows  esat  as  a  function  of  temperature  for  saturation  over  water, 
over  supercooled  water,  and  over  ice.  Instruments  that  measure  the  dew-point  or 
frost-point  temperature  essentially  provide  the  temperature  T  in  (4.1)-(4.4). 
Hence,  Figure  1  is  equivalently  a  plot  of  saturation  vapor  pressure  versus  dew¬ 
point  or  frost-point  temperature,  depending  on  whether  the  surface  in  equilibrium 
is  water  or  ice. 


Figure  1.  Saturation  vapor  pressure  as  a  function  of 
temperature  for  three  regimes:  vapor  over  a  water  surface 
[i.e.,  (4.1)],  vapor  over  supercooled  water  [i.e.,  (4.2)],  and 
vapor  over  ice  [i.e.,  (4.3)].  The  barometric  pressure  is 
assumed  to  be  1000  mb. 
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When  the  water  surface  is  not  pure  but,  for  example,  is  seawater  with  salinity 
S,  the  saturation  vapor  pressure  is  depressed  to  esat(T,S),  where  (e.g.,  Roll  1965, 
p.  262;  List  1984,  p.  373) 

C,*1T,SI  =  1  -  0.000537S  .  (4.5) 

UT) 


Here,  the  salinity  is  in  psu.  Equation  (4.5)  means  that,  for  a  typical  open-ocean 
salinity  of  35  psu,  water  vapor  in  saturation  with  the  surface  has  a  vapor  pressure 
that  is  98.1%  of  the  vapor  pressure  over  a  pure  water  surface  with  the  same 
temperature. 

The  discussion  naturally  turns  now  to  how  to  treat  a  sea  ice  surface.  The 
point  is  often  moot,  however,  because  sea  ice  is  usually  snow  covered,  so  we 
would  just  compute  the  saturation  vapor  pressure  as  that  for  pure  ice  using  (4.3) 
or  (4.4).  If  the  surface  is  truly  bare  sea  ice,  on  the  other  hand,  such  as  for  new  ice 
in  a  freezing  lead  or  summer  sea  ice  when  the  snow  has  all  melted,  we  could  still 
reasonably  just  use  the  saturation  vapor  pressure  relation  for  pure  ice.  Freezing 
seawater  rejects  salt;  consequently,  sea  ice  rarely  has  a  surface  salinity  above 
8  psu.  If  we  assume  that  (4.5)  is  also  appropriate  for  sea  ice — a  reasonable 
assumption — the  depression  in  vapor  pressure  over  sea  ice  with  salinity  8  psu 
would  be  only  about  0.4%.  Most  humidity  sensors  cannot  resolve  such  small 
changes  in  vapor  pressure.  Nevertheless,  if  such  differences  are  important,  (4.5) 
should  be  adequate  for  quantifying  them. 

4.2  Specific  Humidity 

Most  other  water  vapor  variables  are  calculated  from  the  actual  vapor 
pressure  e  and  the  saturation  vapor  pressure  esat.  For  example,  (3.1)  shows  how  to 
compute  the  water  vapor  density  pv  from  the  vapor  pressure.  The  saturation  vapor 
density  pV;Sat  likewise  comes  from  (3.1)  with  esat  substituted  for  e. 

The  specific  humidity  is  defined  as 

Q  =  ~ ,  (4.6) 

Pd  +  Pv 

and  the  saturation  specific  humidity  is 

P  v,sat 


Qsat 


Pd  +  Pv.sat 


(4.7) 
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Notice  that,  from  (3.1)  and  (3.2), 


M„ 


Me 


M,  P 


Ma(P-e)  +  Mwe 


( 


1  - 


1  - 


M„ 


M 


a  J 


0.622(e/P) 

1  -  0.378(e/P)  ' 


(4.8) 


That  is,  the  specific  humidity  also  derives  from  the  vapor  pressure.  In  SI  units, 
specific  humidity  is  given  in  kg  kg~‘ . 


4.3  Mixing  Ratio 

The  mixing  ratio  is  defined  as 


r 


IT 

Pd  ’ 


and  the  saturation  mixing  ratio  is 


(4.9) 


Tat  =  —  •  (4-10) 

Pd 

As  with  (4.8),  we  can  also  write  r  in  terms  of  the  vapor  pressure.  From  (3.1) 
and  (3.2), 


M^e 

Mwe  _  Ma  P  _  0.622(e/P) 

Ma(P-e)“77V“  l-(e/P) 
P 


(4.11) 


4.4  Relative  Humidity 

Officially,  the  relative  humidity  (in  percent)  is  defined  as  the  ratio  of  mixing 
ratios  (Bohren  and  Albrecht  1998,  p.  198;  Glickman  2000,  p.  642f.), 


RH  =100 


4r3 


V  Tat  J 


(4.12) 
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From  the  definitions  of  mixing  ratios,  though,  we  see  that 


RH  =  10(/  6  6sat 


VP  ey 


V  sat  J 


We  can  rearrange  this  as 


RH  =  100 


f  e  ' 


V  ®sat  J 


l-(esa,/P)' 

1  -  (e/P) 


-  100 


1  - 


VSatyL  v 


/  _ 

\ 

- 

esat“ 

e 

1  esa,  e 

l  P 

J 

1  p2 

(4.13) 


(4.14) 


Over  marine  surfaces,  for  example,  esat/P  and  e/P  are  usually  less  than  0.03,  and  e 
is  rarely  less  than  0.5  esat.  Hence,  for  practical  purposes,  we  can  generally  ignore 
the  term  in  square  brackets  in  (4. 14);  (4. 12)  is,  thus,  equivalent  to  the  traditional 
definition  of  relative  humidity, 


RH  =  100 


f  e  ' 


V  esat  J 


(4.15) 


Furthermore,  Bohren  and  Albrecht  (1998,  p.  186)  prefer  (4.15)  to  (4.12)  because 
it  better  reflects  the  physics  of  evaporation  and  condensation  processes.  1 
likewise  use  (4.15)  as  my  definition  of  relative  humidity. 

Table  1  shows  the  relative  humidity  as  a  function  of  air  temperature  and  dew¬ 
point  temperature,  where  I  use  (4.15)  to  define  relative  humidity. 


4.5  Wet-Bulb  Temperature 

The  wet-bulb  temperature,  Twet,  is  another  common  humidity  variable  (e.g., 
Schwerdtfeger  1976,  p  50ff.).  A  wetted  thermometer  will  read  lower  than  the 
ambient  air  temperature  because  of  evaporation  (but  higher  with  condensation), 
and  that  rate  of  evaporation  depends  on  the  mixing  ratio.  At  equilibrium,  a  well- 
ventilated  wet-bulb  thermometer  obeys  (Andreas  1995;  Bohren  and  Albrecht 
1998,  p.  218ff) 


Cp(Ta)(Ta  -Twet)  =  Lv  (Ta  )[rSa,  (Twet  )  “  r]  » 


(4.16) 
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where  cp  =  specific  heat  of  dry  air  at  temperature  Ta 
Lv  =  latent  heat  of  vaporization  at  Ta 
r  =  ambient  mixing  ratio 

rsat  =  saturation  mixing  ratio  at  the  wet-bulb  temperature. 

Rearranging  (4.16)  shows  that  the  wet-bulb  temperature  predicts  the  mixing 
ratio, 


r  =  r  (T  )  -  Cp  (Ta  ^  (T  -  T  ) 

A  Asat  V  wet )  y  /m  \\la  Awet )  ’ 

Lv(ia) 


(4.17) 


Table  1.  Relative  humidity  (RH,  in  percent),  from  (4.15),  as  a  function  of  the  air  temperature  (Ta)  and  the 
dew-point  depression,  where  Td  is  the  dew-point  temperature  (frost-point  temperature  for  Td  less  than 
0°C).  The  barometric  pressure  P  is  assumed  to  be  1000  mb. 
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And  in  light  of  (4. 1 1), 


e 


Aat  (  T  ,,  ) 


f  P  -e  N 


P  -  e 


sat  J 


Ma  (P-e)cp(T.) 

Mw  Lv(Ta)  ‘  wet) 


(4.18) 


which  is  approximately 


®  ®sat('^'wet) 


Pcp(T.) 

0.622  Lv(Ta) 


(Ta-Twet)  . 


(4.19) 


That  is,  the  wet-bulb  temperature  also  predicts  the  vapor  pressure.  Table  2  shows 
the  vapor  pressure  computed  from  (4.19)  for  a  range  of  Ta  and  Ta  -  Twet  values. 

Equations  (4.17)  and  (4.19)  apply  to  perfect  wet-bulb  thermometers. 
Schwerdtfeger  (1976,  p.  5 If.)  and  Brock  and  Richardson  (2001,  p.  94)  describe 
some  second-order  corrections  that  may  be  necessary  to  account  for  ventilation 
rate,  radiative  effects,  and  the  size  of  the  wet  bulb. 
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Table  2.  Vapor  pressure  (e,  in  mb)  from  (4.19)  as  a  function  of  the  air  temperature  (Ta)  and  the  wet-bulb 
depression  for  wet-bulb  temperatures  above  freezing.  The  column  with  Ta  -  Twet  =  0  is  also  the 
saturation  vapor  pressure  (esat)  at  the  indicated  air  temperature;  consequently,  the  ratio  of  e  at 
Ta  -  Twet  >  0  to  esat  is  the  fractional  relative  humidity  (f)  at  Ta  -  Twet.  The  barometric  pressure  P  is 
assumed  to  be  1000  mb. 
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5  WATER  DENSITY 

I  have  found  several  modem  expressions  in  the  literature  for  the  density  of 
pure  water,  pw.  Any  one  would  probably  be  accurate  enough  for  the  purposes  of 
this  handbook.  This  one  is  Gill’s  (1982,  p.  599): 

pw  =  999.842594  +  6.793952xlCT2  T  -  9.095290xl0“3  T2 

■  (5.1) 

+  1.001685X10-4 T3  -  1.120083xl0“6  T4  +  6.536332xl0“9 T5 

It  gives  pw  in  kg  m  for  a  barometric  pressure  of  one  atmosphere  when  T  is 
between  0°  and  40°C. 

To  find  the  density  of  seawater,  psw,  with  temperature  T  and  salinity  S  for 
pressures  near  one  atmosphere,  Gill  (1982,  p.  599)  gives 

psw  =  Pw  +  S(0. 824439  -  4.0899xl(T3  T  +  7.6438xl(T5  T2 
-  8.2467 xlO-7  T3  +  5.3875xl0“9 T4) 

■  (5.2) 

+  S3/2(-5.72466xl0“3  +  1.0227xl0“4 T  -  1.6546xl0“6  T2) 

+  4.8314xl0~4  S2 

This  is  appropriate  for  S  in  [0,  42  psu]  and  T  in  [Tf,  40°C],  where  Tf  is  the 
freezing  point  of  seawater  with  salinity  greater  than  0  psu. 

Gill  (1982,  p.  602)  gives  Tf  as  a  function  of  salinity  for  S  between  0  and 
40  psu;  but  1  prefer  Kester’s  (1974)  formula,  which  is  always  within  0.0 1°C  of 
Gill’s  for  S  in  [1,  40  psu],  because  it  is  easier  to  invert.  Kester  gives  Tf  in  °C  as 

Tf  =  -0.0137  - 5. 1990x10“2S-  7.225x10"5S2  .  (5.3) 

In  turn,  if  we  know  that  the  seawater  is  at  its  freezing  point,  we  can  invert  (5.3)  to 
find  the  corresponding  salinity, 

S  =  -  3.598x  102 

+  6.920x103  [2.7030x10'3  -  2.890xl0‘4  (Tf  +  0.0137)]12  ' 


(5.4) 
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For  supercooled  pure  water  in  the  temperature  interval  [-33°,  0°C], 
Pruppacher  and  Klett  (1997,  p.  87)  give 


pw  =  999.86  +  6.690xl(T2  T  -  8.486xl(T3  T2  +  1.518x1 0-4  T3 
-  6.9984x10“6T4  -  3.6449xl0‘7  T5  -  7.497  xl0“9  T6 


(5.5) 


Here,  again,  pw  is  in  kg  m  3  when  T  is  in  °C.  Pruppacher  and  Klett  have  taken 

(5.5)  directly  from  Hare  and  Sorensen  (1987). 

For  completeness,  I  also  include  in  this  section  an  equation  for  the  density  of 
pure  ice,  pi.  Pruppacher  and  Klett  (1997,  p.  79f.)  give 

p,  =  916.7  -  0.175T  -  5-OxlO-4  T2  ,  (5.6) 

which  yields  pj  in  kg  m  for  T  in  °C.  Pruppacher  and  Klett  claim  that  this  relation 
fits  the  experimental  data  for  T  in  [-180°,  0°C],  On  comparing  the  predictions  of 

(5.6)  with  Hobbs’s  (1974,  p.  348)  tabulation  of  pl5 1  can  further  say  that  (5.6)  is 
accurate  to  better  than  0.5%  over  this  range. 
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6  SPECIFIC  HEAT 

6.1  Specific  Heat  of  Dry  Air 

The  specific  heat  of  air  appears  recurrently  in  studies  of  the  atmospheric 
boundary  layer.  Using  data  from  Hilsenrath  et  al.  (1960),  I  have  obtained  the 
following  polynomial  prediction  for  the  specific  heat  of  dry  air  at  constant 
pressure: 

Cpd  =  1005.60  +  0.01721  IT  +  0.000392T2  .  (6.1) 

This  gives  cpd  in  J  kg~'  °C  1  for  temperature  T  between  -40°  and  40°C  and  for 
barometric  pressures  near  one  atmosphere.  Equation  (6.1)  has  a  minimum  of 
1005.41  J  kg  1  °C at  -21.95°C. 

6.2  Specific  Heat  of  Water  Vapor 

Reid  et  al.  (1987,  p.  656f.,  668)  give  a  polynomial  expression  for  the  specific 
heat  of  water  vapor  at  constant  pressure  for  barometric  pressures  near  one 
atmosphere.  The  temperature  in  their  polynomial  is  in  kelvins,  however,  and  their 
units  of  cpy  are  J  mof1  °C  1  have  therefore  converted  their  polynomial  to 

cpv  =  1858  +  3.820x101T  +  4.220x10“4T2  -1.996x10“7T3  ,  (6.2) 

which  gives  cpv  in  J  kg-1  °C  1  when  T  is  in  °C.  Equation  (6.2)  should  be  accurate 
for  all  near-surface  atmospheric  temperatures. 

6.3  In  the  Turbulent  Sensible  Heat  Flux 

A  frequent  use  for  the  specific  heat  of  air  is  in  finding  the  turbulent  sensible 
heat  flux,  which  is  defined  as 

Hs  -  pacPwt  .  (6.3) 

Elere,  w  is  the  turbulent  fluctuation  in  vertical  wind  velocity,  t  is  the  turbulent 
fluctuation  in  air  temperature,  and  the  overbar  denotes  a  time  average.  The  tildes 
over  the  pa  and  the  cp  terms  in  (6.3)  denote  these  as  general  values  of  the  air 
density  and  specific  heat  of  air  at  constant  pressure  because  confusion  exists  as  to 
which  values  constitute  the  proper  definition  of  sensible  heat  flux. 
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Most  of  the  boundary-layer  community  simply  use  pa,  the  density  of  moist 
air,  for  pa  and  cpd  for  cP .  Businger’s  (1982)  analysis  confirms  that  this  is  proper 
practice  for  micrometeorological  studies.  That  is,  (6.3)  would  be 


Hs  =  Pa  cpd  wt  .  (6.4) 

On  the  other  hand,  when  long  time  intervals  or  global  averages  define  the 
scope  of  the  study — for  example,  in  balancing  the  hydrological  cycle — Businger 
(1982)  explains  that  the  reference  temperature  for  enthalpy  transfer  must  be 
chosen  carefully.  When  such  “careful  bookkeeping  of  the  energy”  is  necessary, 
Businger  concludes  that  the  sensible  heat  flux  must  be  expressed  as 

Hs  =  (PdcPd  +Pvcpv)wt  •  (6-5) 

Fuehrer  and  Friehe’s  [2002,  Eq.  (1 14)]  extensive  thermodynamic  analysis 
yields  essentially  this  same  result,  although  they  give  it  as  general  result — not 
one  necessary  only  for  large  areal  averages.  They  also  include  two  other  small 
terms  in  (6.5)  that  are  associated  with  the  flux  of  water  vapor  between 
temperature  reference  states. 

As  a  practical  exercise,  1  can  check  how  different  (6.4)  and  (6.5)  are.  We  can 
convert  (6.5)  into  a  form  similar  to  (6.4), 


wt  .  (6.6) 

Using  (6.1)  and  (6.2),  we  see  that  (cpv  -  cpd)/cpd  ranges  from  0.861  to  0.833  for  air 
temperatures  between  -40°  and  40°C.  Flence,  1  approximate  (6.6)  as  (cf.  Larsen 
and  Busch  1974) 


Hs  =  pac 


a  pd 


1  +  Q 


rc 

pv  pd 


^pd  J 


Hs  -  Pa cpd (1  +  0.85Q)wt  .  (6.7) 

As  1  mentioned  earlier,  Q  is  seldom  larger  than  0.035  kg  kg~’  in  the  natural 
atmosphere,  and  it  would  attain  this  value,  probably,  only  over  the  tropical  ocean. 
Flence,  (6.5)  yields  values  of  the  sensible  heat  flux  that  are,  at  most,  3%  larger  than 
the  values  from  (6.4).  In  many  applications,  this  is  a  negligible  difference,  since 
wt  usually  is  measured  to  no  better  than  ±10%.  And  at  temperatures  below  0°C, 
where  Q  is  less  than  0.004  kg  kg-1,  (6.4)  and  (6.5)  differ  by  much  less  than  1%. 
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6.4  Specific  Heat  of  Water 

Both  temperature  and  salinity  influence  the  specific  heat  of  seawater  at 
constant  pressure,  cps.  Neumann  and  Pierson  (1966,  p.  47)  tabulate  values  of  cps 
that  come  from  data  collected  by  Cox  and  Smith  (1959).  Home  (1969,  p.  68) 
gives  a  functional  expression  for  cps  in  terms  of  temperature  and  salinity  that  he 
converted  from  a  similar  expression  that  Bromley  et  al.  (1967)  deduced  from 
their  own  measurements. 

Millero  et  al.  (1973)  also  reported  measurements  of  cps  and  fitted  an  equation 
to  these  measurements  in  terms  of  temperature  and  chlorinity.  Gill  (1982,  p.  601) 
converted  this  to  an  expression  in  terms  of  temperature  and  salinity.  The 
following  is  essentially  Gill’s  equation,  although  I  have  modified  it  slightly  to 
better  represent  the  precision  that  Millero  et  al.  implied  in  their  fitting 
coefficients: 


Cps  =  Cpw  +  S(-7.6444  +  0.10728T  -  1.384xl0^3  T2) 
+  S3/2  (0.177  -  4.08 xlO-3  T  +  5.35xl0"s  T2) 

Here,  cpw  is  the  specific  heat  of  pure  water, 


(6.8) 


Cpw  =  4217.4  -  3.720283T  +  0.1412855T2 

-  2.654387x10'3T3  +  2.093236xl0~5  T4  ' 


(6.9) 


In  (6.8)  and  (6.9),  T  is  in  °C  and  ranges  from  Tf  to  40°C;  in  (6.8),  S  is  in  psu  and 
ranges  from  0  to  40  psu.  The  pressure  is  assumed  to  be  one  atmosphere.  Table  3 
shows  values  of  cps  that  result  from  (6.8). 

The  values  in  Table  3  are  quite  compatible  with  the  values  in  Neumann  and 
Pierson’s  (1966,  p.  47)  table,  which  are  based  on  Cox  and  Smith’s  (1959)  data. 
The  values  in  Table  3,  however,  have  some  systematic  differences  from  the 
results  from  Bromley  et  al.  (1967;  also  Home  1969,  p.  67£),  especially  for  fresh 
water  in  the  lower  part  of  the  temperature  range.  For  example,  Bromley  et  al. 
report  cpw  =  4207  J  kg-1  °C  for  T  =  0°C,  while  Table  3  gives  4217  J  kg-1  °C. 
Therefore,  1  conclude  that  (6.8)  and  (6.9)  are  accurate  to  roughly  ±3  J  kg-1  °C  (cf. 
Bromley  et  al.  1967;  Millero  et  al.  1973). 
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Table  3.  Specific  heat  of  seawater  at  constant  pressure,  cps  (in 
J  kg1  °C_1),  as  a  function  of  temperature  and  salinity.  These  results 
come  from  (6.8).  The  barometric  pressure  is  assumed  to  be  one 
atmosphere. 
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6.5  Specific  Heat  of  Ice 

Murphy  and  Koop  (2005)  give  a  new  expression  for  the  specific  heat  of  ice  for 
temperatures  down  to  20  K.  Because  their  units,  however,  are  J  mol  1  K  I  convert 
their  expression  to  predict  cpi  in  J  kg-1  K  1  using  Mw  =  1 8.0 1 5  x  1 0  3  kg  mol1.  The 
result  is 


cpi  =  -114.19  + 


L1288T  +  3.421T  exp 


-(T/125.1)2 


(6.10) 


in  which  T  must  be  in  ke Ivins.  Murphy  and  Koop  use  (6. 10),  for  example,  to 
derive  their  expression  for  the  saturation  vapor  pressure  over  ice,  (4.4)  above. 

Figure  2  compares  values  for  the  specific  heats  of  air,  pure  water,  ice,  and 
water  vapor  for  temperatures  between  -40°  and  40°C. 
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Temperature  (°C) 

Figure  2.  Specific  heats  at  constant  pressure  of  dry  air 
[from  (6.1)],  water  vapor  [from  (6.2)],  pure  water  [from 
(6.9)],  and  pure  ice  [from  (6.10)].  Notice,  because  the 
vertical  axis  has  such  a  large  range,  the  small  variations 
in  cpd,  cpv,  and  cpw  are  not  obvious. 
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7  LATENT  HEAT 

Fleagle  and  Businger  (1980,  p.  113)  give  the  following  equations  for  the 
latent  heats  associated  with  the  phase  transitions  of  water  molecules: 

Latent  heat  of  vaporization, 

Lv  =  (25.00  -  0.02274T)xl05  ;  (7.1) 

Latent  heat  of  fusion, 

Lf  =  3.34x10s  ;  (7.2) 

Latent  heat  of  sublimation, 

Ls  =  (28.34  -  0.00149T)xl05  .  (7.3) 

Each  of  these  gives  the  latent  heat  in  J  kg  1  when  the  temperature  T  is  in  °C. 

The  value  for  Lf  in  (7.2)  is  appropriate  only  near  0°C  because  Lf  decreases 
with  decreasing  temperatures.  Hobbs’s  (1974,  p.  361)  tabulation  of  Lf  for 
temperatures  below  0°C  disagrees  dramatically  with  a  similar  tabulation  in  the 
Smithsonian  Meteorological  Tables  (List  1984,  p.  343),  however.  The  functions 
for  Lv  and  Ls,  on  the  other  hand,  agree  quite  well  with  the  Smithsonian 
tabulation.  The  Lv  values  predicted  by  (7.1)  are  within  0.3%  of  the  Smithsonian 
values  for  temperatures  from  0°  to  60°C,  and  the  Ls  values  from  (7.3)  are  within 
0.2%  of  the  Smithsonian  values  for  temperatures  between  -50°  and  0°C. 
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8  SURFACE  TENSION  OF  WATER 

Vargaftik  et  al.  (1983)  tabulate  consensus  values  for  the  surface  tension  of 
pure  water  for  temperatures  between  0°C  and  the  critical  temperature,  374°C 
(Bohren  and  Albrecht  1998,  p.  207ff.).  They  also  develop  from  these  data  the 
following  relation  for  computing  the  surface  tension  of  pure  water: 


a 


w 


f 

0.2358 

V 


374.00-T 

647.15 


n  1.256 


f 

1  -  0.625 

V 


374.00-T" 
647.15  , 


(8.1) 


This  gives  aw  in  J  m  2  when  the  water  temperature  T  is  in  °C.  This  equation 
predicts  exactly  the  values  of  surface  tension  that  Lide  (2001,  p.  6-3)  tabulates 
and  matches  the  values  in  Batchelor  (1970,  p.  597)  to  within  about  0.1%. 

For  temperatures  between  -45°  and  0°C,  Pruppacher  and  Klett  (1997,  p.  130) 
suggest  computing  the  surface  tension  of  pure  water  from 


gw  =  7.593xl0~2  +  1.15X10”4  T  +  6.818xlO'5T2  +  6.511xlO_6T3 

,  (8.2) 

+  2.933xl0“7  T4  +  6.283xl0“9  T5  +  5.285xl0“n  T6 

which  again  gives  gw  in  J  m  2  for  T  in  °C.  Unfortunately,  (8.1)  and  (8.2)  do  not 
meet  at  0°C.  Equation  (8.1)  predicts  7.656  x  10  2  J  nT2  there,  while  (8.2)  gives 
7.593  x  10  2  J  nT2.  The  former  value  is  probably  the  more  accurate  one. 

Pruppacher  and  Klett  (1978,  p.  107)  likewise  give  an  expression  for  the 
surface  tension  of  an  interface  between  water  vapor  and  saline  water.  That 
expression,  however,  quantifies  the  salinity  in  terms  of  ms/mw,  where  mw  is  the 
mass  of  pure  water  per  unit  volume  and  ms  is  the  mass  of  dissolved  salt  in  the 
volume.  Since  the  definition  of  salinity  is 


s 


m„ 


m,„  +  m„ 


(8.3) 


we  see  that 


JEi 

m.„ 


1  —  s  ’ 


(8.4) 
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where  s  is  the  fractional  salinity  [s  =  S(in  psu)/1000]. 

Hence,  I  convert  Pruppacher  and  Klett’s  (1978)  expression  for  the  surface 
tension  at  a  seawater  interface  to 


ci, 


=  o,„  +  2.77x10“ 


1  —  s 


(8.5) 


where  aw  comes  from  (8.1)  or  (8.2),  depending  on  the  water  temperature.  The 
coefficient  multiplying  the  salinity  term  in  (8.5)  is  virtually  the  same  value  that 
Hanel  (1976)  recommends.  Equation  (8.5)  should  be  accurate  for  temperatures  in 
[-25°,  40°C]  and  S  in  [0,  260  psu]. 
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9  MEAN  FREE  PATHS  OF  AIR  AND  WATER  VAPOR 
MOLECULES 

The  mean  free  path  of  a  gas  molecule  is  an  estimate  of  the  distance  it  travels 
between  collisions  with  other  molecules  in  the  gas.  Starting  with  equations  that 
Wagner  [1982,  Eqs.  (5.55),  (5.56)]  gives,  I  derive  these  expressions  for  the  mean 
free  paths  of  air  molecules  (Xa)  and  water  vapor  molecules  (/.v)  in  air: 


X 


a 


kT 

V2n(P-e)d;  ’ 


(9.1) 


/L622  Tt(P-e)(da  +  dw)2 


where  k  =  Boltzmann  constant 
T  =  temperature  (K) 

P  =  barometric  pressure  (Pa) 
e  =  water  vapor  pressure  (Pa) 
da  =  diameter  of  an  air  molecule 
dw  =  diameter  of  a  water  molecule. 

Equation  (9.1)  is  essentially  the  same  as  Reif  s  (1967)  expression  for  the  mean 
free  path  of  air  molecules, 


X 


a 


kT 

V2  xtPda 


(9.3) 


In  (9.1)  and  (9.2),  Xa  and  Xv  quantify  the  distance  between  consecutive 
collisions,  not  the  distance  between  collisions  between  molecules  of  the  same 
species.  That  is,  a  water  vapor  molecule  will  likely  hit  an  air  molecule  next;  Xv 
estimates  the  distance  from  the  previous  to  this  next  collision. 

Using  values  for  da  and  dw  of  3.7  x  1CT10  m  and  1.65  x  1CT10  m,  respectively*, 
and  expressing  P  and  e  in  millibars,  I  simplify  (9.1)  and  (9.2)  to 


Personal  communication,  James  H.  Cragin,  CRREL,  2000. 
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.  2.3xl(T7  T 

D 

P-e 

(9.4) 

.  4.8x1  O'7  T 

=  P-e  ' 

(9.5) 

In  these,  A.a  and  /-v  are  in  meters  when  T  is  in  kelvins  and  P  and  e 

are  in  millibars. 

For  example,  when  T  =  293  K  and  P  -  e  =  1000  mb. 

A.a  =  6.7 xl0~8  m  , 

(9.6) 

Xv  -  1.4xl0“7  m  . 

(9.7) 

For  comparison,  Pruppacher  and  Klett  (1978,  p.  323)  state  without  proof  that 
=  6.6  x  1CT8  m  for  P  =  P0  and  T  =  293.15  K.  Bohren  and  Albrecht  (1998,  p. 
68)  estimate  A,a  =  1 0  7  m  for  pressures  near  one  atmosphere. 
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10  MOLECULAR  VISCOSITY 

10.1  Air 

Hilsenrath  et  al.  (1960,  p.  10)  give  the  following  expression  for  the  dynamic 
viscosity  of  air: 


1.458xlQ~6  T3/2 
T  +  110.4 


(10.1) 


Here,  r\d  is  in  kg  m  1  s  1  when  T  is  in  kelvins,  and  they  suggest  that  this  relation 
is  accurate  for  a  pressure  of  one  atmosphere  for  temperatures  between  1 00  and 
1900  K. 

The  kinematic  viscosity  of  air,  va,  however,  occurs  more  commonly  in 
boundary-layer  studies.  We  can  obtain  this  from  pa  as 


va  =  ^  •  (10.2) 

Pd 

1  have  combined  this  definition,  (2.2),  and  (10.1)  to  obtain  a  polynomial 
expression  for  predicting  va  for  the  temperature  range  [-50°,  50°C]; 

va  =  1.326xl0“5(l  +  6.542x10^3T  +  8.301x10“6T2  -  4.840x10“9T3).  (10.3) 

Here,  va  is  in  m2  s~’,  and  T  is  in  °C. 

Figure  3  shows  (10.3)  plotted  as  a  function  of  temperature;  Table  4  lists  these 
values.  The  predictions  from  (10.3)  agree  to  three  significant  figures  with  values 
tabulated  in  Goldstein  (1965,  p.  7)  and  Batchelor  (1970,  p.  594). 


10.2  Water 

Reid  et  al.  (1987,  p.  441,  455)  give  the  following  expression  for  the  dynamic 
viscosity  of  pure  water: 

r|w  =  10~3  exp (-24.71  +  4.209xl03  T'1  +  4.527xl0“2  T  -  3.376xl0“5  T2) .  (10.4) 
As  with  (10. 1),  here  pw  is  in  kg  m  1  s  ',  and  the  temperature  is  in  kelvins. 
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Temperature  (°C) 

Figure  3.  Molecular  values  of  the  kinematic  viscosity  (va)  and 
thermal  diffusivity  (D)  of  air  and  of  the  water  vapor  diffusivity 
(Dv)  in  air  at  a  barometric  pressure  of  1000  mb.  The  calcula¬ 
tions  of  viscosity  and  thermal  diffusivity  assume  dry  air. 

Because  Neumann  and  Pierson  (1966,  p.  52)  suggest  that  salinity  and 
pressure  affect  the  dynamic  viscosity  of  water  only  slightly,  1  compute  the 
kinematic  viscosity  of  pure  water  (vw)  and  seawater  (vsw)  from 

vw=^  (10.5) 

Pw 

and 

vsw  =  ^  •  (10.6) 

Psw 

In  these,  vw  and  vsw  are  in  m2  s-1;  and  pw  and  psw  come  from  (5.1)  and  (5.2), 
respectively. 

As  examples  of  these  calculations,  for  T  =  20°C,  S  =  34  psu,  and  pressures 
near  one  atmosphere,  I  find  pw  =  1.018  x  1 0  kg  m  1  s  ',  pw  =  998.2  kg  m  \ 
psw  =  1024.0  kg  mf3,  vw  =  1.019  x  1(T6  m2  s“\  and  vsw  =  0.994  x  10‘6  m2  s-1. 

In  particular,  these  pw  and  vw  values  and  values  1  have  calculated  at  other 
temperatures  are  within  3%  of  values  tabulated  in  Home  (1969,  p.  486)  and 
Batchelor  (1970,  p.  597). 
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Table  4.  Molecular  values  of  the  kinematic  viscosity  (va) 
and  thermal  diffusivity  (D)  of  air,  the  diffusivity  of  water 
vapor  in  air  (Dv),  and  the  Prandtl  (Pr)  and  Schmidt  (Sc) 
numbers.  Calculations  of  the  thermal  diffusivity  assume 
dry  air.  The  barometric  pressure  is  assumed  to  be 
1000  mb. 
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1 1  THERMAL  CONDUCTIVITY  AND 
THERMAL  DIFFUSIVITY  OF  AIR 

Hilsenrath  et  al.  (1960,  p.  70)  also  tabulate  the  thermal  conductivity  of  air,  ka. 
1  have  fitted  their  data  for  temperatures  between  -193°  and  277°C  with  the 
following  polynomial: 

ka  =  2.4 1  lxl 0-2 (l  +  3.309x10~3T  -  1.441xlO“6T2)  .  (11.1) 

This  gives  ka  in  W  m  1  °C  1  for  temperature  T  in  °C. 

The  thermal  diffusivity  D  is  analogous  as  a  molecular  transport  variable  to 
the  kinematic  viscosity  va.  We  define  D  as 

D  =  — ,  (11.2) 

Pd  cp 

which  has  units  of  m2  s  Figure  3  plots  D  as  a  function  of  temperature,  while 
Table  4  lists  values  for  normal  atmospheric  boundary  layer  temperatures. 

Notice  that,  in  general,  pa  and  cp  in  (11.2)  should  include  the  effects  of 
atmospheric  water  vapor.  1  have  ignored  those  small  effects  in  creating  Figure  3 
and  Table  4  because  (11.1)  predicts  the  thermal  conductivity  of  dry  air.  1  presume 
that  including  water  vapor  effects  in  ka  would  tend  to  offset  the  effects  of  water 
vapor  in  pa  and  cp  and  would,  thus,  yield  a  thermal  conductivity  D  comparable  to 
the  dry-air  value  that  1  have  calculated. 


Physical  Constants  and  Functions  for  Atmospheric  Boundary  Layer  Studies 


29 


12  MOLECULAR  DIFFUSIVITIES  OF  GASES  IN  AIR 


12.1  Water  Vapor 

Hall  and  Pruppacher  (1976)  developed  the  following  expression  for  the 
molecular  diffusivity  of  water  vapor  in  air: 


D„  =  2.11x10“ 


7  T  V-94 
vToy 


f  D  A 


V  r  j 


(12.1) 


This  gives  Dv  in  m2  s  1  for  temperature  T  in  kelvins  and  pressure  P  in  millibars. 
Hall  and  Pruppacher  explain  that  no  good  measurements  of  Dv  exist  for 
temperatures  below  0°C.  Therefore,  they  obtain  (12.1)  by  extrapolating 
measurements  above  freezing  to  subfreezing  temperatures.  Still,  Hall  and 
Pruppacher  claim  that  (12.1)  applies  over  the  temperature  interval  [-80°,  40°C], 
Pruppacher  and  Klett  (1978,  p.  413)  originally  recommended  (12.1),  and 
Pruppacher  and  Klett  (1997,  p.  503)  still  do.  Figure  3  and  Table  4  compare 
values  of  Dv  with  va  and  D. 

The  Prandtl  (Pr)  and  Schmidt  (Sc)  numbers  compare  the  relative  efficiencies 
of  molecular  exchange  processes.  The  molecular  Prandtl  number  for  air  is 


Pr  =  — 

D 


and  the  molecular  Schmidt  number  for  air  is 


(12.2) 


Sc  =  —  .  (12.3) 

Dv 

We  can  compute  these  from  (10.3),  (11.2),  and  (12.1).  Table  4  also  lists  values  of 
Pr  and  Sc  for  P  =  1000  mb  and  temperatures  in  the  range  [-40°,  40°C], 

12.2  Other  Atmospheric  Gases 

The  molecular  diffusivities  of  atmospheric  trace  gases  such  as  carbon 
dioxide,  methane,  and  nitrous  oxide  also  occur  in  atmospheric  boundary  layer 
research.  Reid  et  al.  (1987,  p.  587)  give  a  semi-empirical  expression,  which  they 
adapted  from  Fuller  et  al.  (1969),  to  predict  how  the  diffusivities  of  these  trace 
atmospheric  gases  depend  on  temperature  and  pressure. 
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Table  5  lists  diffusivities  in  air  for  several  environmentally  important  gases.  I 
have  computed  some  of  these  (labeled  “Dg0,  Reid  et  al”)  for  T0  =  273.15  K  and 
Po  =  1013.25  mb  from  equation  (1 1-4.4)  and  Table  11.1  in  Reid  et  al.  (1987).  For 
comparison,  1  also  include  in  Table  5  diffusivity  values  (labeled  DgT)  at 
temperatures  of  0°C  (273.15  K)  or  25°C  (298.15  K)  that  Thibodeaux  (1979, 

1996)  tabulates.  Thibodeaux,  however,  does  not  mention  the  barometric  pressure 
corresponding  to  his  values;  I  therefore  assume  it  is  approximately  one 
atmosphere. 

Equation  (1 1-4.4)  in  Reid  et  al.  (1987)  implies  a  general  expression  for  how 
the  diffusivities  of  gases  in  air  depend  on  temperature  and  pressure, 


Ds  =  DS0 


f  T  Y-75  ( 


vToy 


v  r  J 


(12.4) 


In  this,  Dg  is  in  m2  s_1,  T  is  in  kelvins,  P  is  in  millibars,  and  Dg0  is  the  value  at  T0 
and  P0  in  Table  5.  Alternatively,  we  could  just  as  well  substitute  Thibodeaux’s 
(1979,  1996)  values  for  DgT  in  place  of  Dg0  in  (12.4),  replace  T0  with  the 
tabulated  temperature,  and  assume  that  Thibodeaux’s  listed  values  all  correspond 
to  pressure  P0. 

Massman  (1998)  also  reviews  the  diffusivities  of  various  gases  in  air  in  the 
context  of  an  equation  like  (12.4)  but  assumes  that  the  exponent  is  1.81: 


Dg  =  °8o 


4  T  Y'81  f 


VTo  J 


V  r  J 


(12.5) 


Table  5  lists  the  values  that  he  recommends  for  Dg0.  Coincidently,  the 
Smithsonian  Meteorological  Tables  (List  1984,  p.  395)  also  suggests  a  relation 
like  (12.5)  for  the  diffusivities  of  gases  in  air. 

Comparing  the  DgT  values  reported  at  273.15  K  directly  with  the  Dg0  values 
in  Table  5  or  using  (12.4)  or  (12.5)  to  convert  the  Reid  et  al.  (1987),  Massman 
(1998),  or  Thibodeaux  (1979,  1996)  values  to  the  same  temperature  lets  us 
evaluate  the  uncertainty  in  these  gas  diffusivities. 

For  example,  the  four  diffusivities  listed  in  Table  5  for  carbon  dioxide  have  a 
spread  of  only  0.06  x  ]  0  nr  s  1  at  0°C.  The  diffusivities  for  ammonia,  on  the 
other  hand,  have  a  spread  of  about  0.5  x  ]  ()  m2  s  1  at  0°C. 

Table  5  also  includes  four  new  estimates  of  the  water  vapor  diffusivity  to 
compare  with  (12.1).  The  two  values  from  Thibodeaux  (1979,  1996)  both  predict 
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Table  5.  Molecular  diffusivities  of  various  gases  in  air.  The  columns  labeled 
Dg0  show  predictions  based  on  equation  (11-4.4)  and  Table  11.1  in  Reid  et 
al.  (1987)  and  recommendations  by  Massman  (1998).  These  values  are 
appropriate  at  temperature  T0  =  273.15  K  and  pressure  P0  =  1013.25  mb.  The 
DgT  column  shows  values  tabulated  in  Thibodeaux  (1979,  Table  C.8, 
superscript  1)  or  Thibodeaux  (1996,  Table  C.6,  superscript  2)  for  the 
temperatures  indicated. 


Gas 

105  Dg0 

(m2  s -1) 

Temperature 

(K) 

105  DgT 
(mV1) 

Reid  et  al. 

Massman 

Ammonia,  NH3 

1.88 

1.98 

273.15 

2.161 

298.15 

2.82 

Bromine,  Br2 

0.971 

298.15 

1.002 

Carbon  dioxide,  C02 

1.35 

1.38 

273.15 

1.381 

298.15 

1.642 

Carbon  monoxide,  CO 

1.71 

1.81 

298.15 

2.032 

Hydrogen  sulfide,  HS 

298.15 

1.662 

Methane,  CH4 

1.95 

273.15 

1.62 

Nitrogen  dioxide,  N02 

1.36 

Nitric  oxide,  NO 

1.80 

298.15 

2.042 

Nitrous  oxide,  N20 

1.22 

1.44 

298.15 

1.552 

Ozone,  03 

1.44 

Sulfur  hexafluoride,  SF6 

0.795 

Sulfur  dioxide,  S02 

1.08 

1.09 

273.15 

1.032 

Water  vapor,  H20 

2.15 

2.18 

273.15 

2.201 

298.15 

2.562 

Dv  =  2.20  x  10  5  m2  s  1  at  0°C.  The  comparable  value  that  1  predict  from  Reid  et 
al.  (1987)  is  Dv  =  2.15  x  10  5  m2  s_1,  while  Massman  (1998)  recommends 
Dv  =  2.18  x  10~5  m2  s-1.  Hall  and  Pruppacher’s  (1976)  equation,  (12.1),  gives 
Dv  =  2. 1 1  x  10~5  m2  s-1  at  0°C  and  1013.25  mb.  Consequently,  1  conclude  that  the 
diffusivities  for  the  gases  listed  in  Table  5  are  typically  known  to  an  accuracy  of 
±3%.  Massman  similarly  concludes  that  most  of  the  diffusivities  listed  in  Table  5 
have  an  absolute  uncertainty  that  is  no  more  than  5-9%,  though  the  tabulated 
diffusivities  for  03,  NO,  and  N03  might  be  uncertain  by  25%  because  of  the 
paucity  of  data  for  these  molecules. 
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13  EFFECTS  OF  SURFACE  CURVATURE  ON  ka  AND  Dv 

The  transfers  of  heat  and  moisture  at  the  surface  of  small  atmospheric 
particles,  like  cloud  droplets,  sea  spray  droplets,  snowflakes,  and  aerosols,  cannot 
strictly  be  parameterized  in  terms  of  the  ka  and  Dv  values  given  above.  Because 
of  the  extreme  curvature  of  the  surface  of  these  small  particles,  the  air  around 
them  no  longer  behaves  as  a  continuum.  This  is  the  so-called  Kelvin  effect  (e.g., 
Pruppacher  and  Klett  1997,  p.  170).  Likewise,  surface  curvature  will  also  affect 
the  transfers  of  heat  and  moisture  around  the  snow  grains  in  a  snowpack,  which 
is  always  porous. 

Pruppacher  and  Klett  (1997)  present  equations  to  account  for  how  these  cur¬ 
vature  effects  modify  ka  and  Dv  (also  Andreas  1989,  1990,  1995).  For  predicting 
the  effects  of  curvature  on  ka,  they  give  (Pruppacher  and  Klett  1997,  p.  509) 


k'  = 


+ 


/  \l/2 

'  27tM„ 


r0+AT  r0aTpacp^  RT  ) 


(13.1) 


where  ka 

r0 

Ay 

aT 

T 


thermal  conductivity  modified  for  curvature  effects  (W  m  1  K  ') 

particle  radius  (m) 

empirical  constant  (=2.16  x  ]  ()  7  m) 

empirical  constant  (=  0.7) 

temperature  (K). 


In  (13.1),  we  recognize  ka/pacp  from  (11.2)  as  the  thermal  diffusivity  D. 
Hence,  if  we  divide  (13.1)  by  pacp,  we  get  an  analogous  expression  for  curvature 
effects  on  the  thermal  diffusivity, 


D' 


D 


r0  +  At 


+ 


D 


/  \l/2 

^2nM^ 


(13.2) 


This  predicts  the  modified  thermal  diffusivity  in  m2  s  1 . 

Similarly,  Pruppacher  and  Klett  (1997,  p.  506)  give  the  following  equation 
for  predicting  how  surface  curvature  influences  the  vapor  diffusivity  around 
small  particles: 
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d:  = 


D 


(13.3) 


where  □'  =  water  vapor  diffusivity  modified  for  curvature  effects  (m2  s  ') 
ac  =  empirical  constant  (=  0.036) 

Av  =  empirical  constant  (=  1.3A,a  =  8.7  x  10~8  m  for  typical  sea-level 
temperature  and  pressure). 

Figure  4  compares  k(  with  ka  and  D'  with  Dv  for  particles  with  radii 
between  0.1  and  1000  pm,  which  is  a  typical  size  range  for  sea  spray  droplets  (cf. 
Andreas  1989).  According  to  this  figure,  curvature  effects  significantly  influence 
sensible  heat  transfer  only  around  particles  with  radii  less  than  about  5  pm.  In 
contrast,  surface  curvature  significantly  decreases  the  rate  of  vapor  diffusion 
around  particles  with  radii  up  to  about  200  pm. 


Figure  4.  Effects  of  surface  curvature  on  the  thermal 
conductivity  of  air  [k'a,  from  (13.1)]  and  on  the  water 
vapor  diffusivity  in  air  [D'v,  from  (13.3)].  The  plot  also 
shows  for  comparison  the  unmodified  values,  ka  [from 
(11.1)]  and  Dv  [from  (12.1)].  The  air  temperature  is 
assumed  to  be  20°C,  and  the  barometric  pressure  is 
1000  mb. 
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14  WIND  FORCE  SCALES 

14.1  Beaufort  Scale 

A  common  way  to  describe  wind  and  sea  state  is  with  the  Beaufort  Scale.  In 
the  early  nineteenth  century,  Admiral  Sir  Francis  Beaufort  developed  a  scale  for 
wind  force  based  on  the  behavior  of  sailing  ships.  The  British  Admiralty  later 
adopted  this  scale,  and  Bowditch  (1977,  p.  1059)  helped  popularize  it. 

In  effect,  the  Beaufort  number  or  force  B  is  related  to  the  wind  speed  (in 
m  s  ')  at  a  standard  reference  height  of  10  m,  Uio,  through  (List  1984,  p.  119; 
Strangeways  200 1 ) 

U10  =  0.836B372  .  (14.1) 

But  the  key  feature  of  the  Beaufort  Scale  is  that  it  associates  Ui0  and  B  with  a 
description  of  wind  and  sea  state  and,  thus,  provides  an  estimate  of  Uio  from 
visual  observations  alone. 

Subsequently,  the  Beaufort  Scale  was  adapted  to  use  over  land  (e.g.,  Bowditch 
1977,  p.  1059;  List  1984,  p.  1 19).  Table  6  shows  the  Beaufort  Scale  and  includes 
descriptions  of  conditions  for  a  given  Beaufort  force  over  both  land  and  sea. 

14.2  Saffir-Simpson  Scale 

The  Beaufort  Scale  classifies  all  ocean  storms  with  surface-level  winds  above 
32.7  m  s  1  as  hurricanes.  But  conditions  at  sea  and  when  the  storm  comes  ashore 
vary  widely  depending  on  the  wind  speed.  The  Saffir-Simpson  Scale,  developed 
by  Herbert  Saffir  and  Bob  Simpson,  further  divides  hurricanes  into  five  categories. 
Table  7  shows  the  Saffir-Simpson  Scale. 

At  sea,  storms  are  assigned  to  a  Saffir-Simpson  category  on  the  basis  of  their 
maximum  surface-level  wind  speed  and  minimum  central  pressure.  Storms  may 
change  category  during  their  lifetime  as  they  intensify  or  degrade.  The  “Storm 
Surge”  listed  in  Table  7  is  the  height  of  the  ocean  wave  that  comes  ashore  ahead 
of  the  storm.  The  values  shown  give  a  typical  range;  the  actual  storm  surge  will 
depend  on  the  slope  of  the  continental  shelf. 

The  main  relevance  of  the  Saffir-Simpson  Scale  is  that  it  attempts  to  forecast 
flooding  and  damage  if  the  storm  does  move  onshore.  The  “Effects”  column  in 
the  table  lists  these  predictions.  Effects  range  from  minor  for  a  Category  1  storm 
to  catastrophic,  as  they  were  with  Mitch,  a  1998  Category  5  hurricane  that  killed 
over  800  people  in  Honduras  and  Nicaragua. 


Table  6.  Beaufort  Scale,  with  the  associated  wind  speed  ranges  for  each  Beaufort  force  in  meters  per  second,  knots,  and  miles  per 
hour.  H1/3  is  the  significant  wave  height,  the  average  height  of  the  highest  one-third  of  all  waves  occurring  during  a  period  (Kinsman 
1965,  p.  302f,  390f.). 


Force 

Wind 

description 

Uio 

H-l/3 

(m) 

Over  the  sea 

Over  land 

(ms  ) 

(knots) 

(mph) 

0 

Calm 

0.0-0. 2 

<1 

<1 

0 

Sea  like  a  mirror 

Calm;  smoke  rises  vertically 

1 

Light  air 

0.3-1 .5 

1-3 

1-3 

CNJ 

o 

1 

o 

Ripples  with  appearance  of  scales;  no  foam  crests 

Smoke  drift  indicates  wind  direction; 
vanes  do  not  move 

2 

Light 

breeze 

1. 6-3.3 

4-6 

4-7 

0.3-0. 5 

Small  wavelets;  crests  have  glassy  appearance  but 
do  not  break 

Wind  felt  on  face;  leaves  rustle;  vanes 
begin  to  move 

3 

Gentle 

breeze 

3. 4-5. 4 

7-10 

8-12 

0.6-1. 0 

Large  wavelets;  crests  begin  to  break;  scattered 
whitecaps 

Leaves  and  twigs  in  constant  motion; 
light  flags  extended 

4 

Moderate 

breeze 

5. 5-7. 9 

11-16 

13-18 

1.5 

Small  waves  becoming  longer;  numerous  whitecaps 

Dust,  leaves,  and  loose  paper  raised; 
small  branches  move;  flags  flap 

5 

Fresh 

breeze 

8.0-10.7 

17-21 

19-24 

2.0 

Moderate  waves  taking  longer  form;  many  whitecaps 
and  chance  of  some  spray 

Small  trees  in  leaf  begin  to  sway; 
whitecaps  on  inland  waters 

6 

Strong 

breeze 

10.8-13.8 

22-27 

25-31 

3.5 

Large  waves  forming;  white  foam  crests  extensive, 
and  spray  probable 

Larger  branches  of  trees  in  motion; 
flags  pop;  whistling  in  wires;  umbrellas 
unstable 

7 

Moderate 

gale 

13.9-17.1 

28-33 

32-38 

5.0 

Sea  heaps  up,  and  white  foam  from  breaking  waves 
begins  to  be  blown  in  streaks;  spindrift  appears 

Whole  trees  in  motion;  resistance  felt 
in  walking  against  the  wind 

8 

Fresh  gale 

17.2-20.7 

34-40 

39-46 

7.5 

Moderately  high  waves  of  greater  length;  edges  of 
crests  break  into  spindrift;  foam  is  blown  in  well- 
marked  streaks 

Twigs  and  small  branches  broken; 
progress  generally  impeded 

9 

Strong  gale 

20.8-24.4 

41-47 

47-54 

9.5 

High  waves;  dense  streaks  of  foam;  sea  begins  to 
roll;  spray  may  reduce  visibility 

Slight  structural  damage  occurs;  slate 
blown  from  roofs 

10 

Whole  gale 

24.5-28.4 

48-55 

55-63 

12 

Very  high  waves  with  overhanging  crests;  sea 
surface  takes  on  white  appearance  as  foam  in  great 
patches  is  blown  in  very  dense  streaks;  rolling  sea  is 
heavy;  visibility  reduced 

Seldom  experienced  on  land;  trees 
broken  or  uprooted;  considerable 
structural  damage  occurs 

11 

Storm 

28.5-32.6 

56-63 

64-72 

15 

Exceptionally  high  waves;  sea  covered  with  long 
white  patches  of  foam;  small  and  medium-sized 
ships  might  be  lost  to  view  behind  waves;  visibility 
further  reduced 

Very  rarely  experienced  on  land; 
usually  accompanied  by  widespread 
damage 

12 

Hurricane 

>32.7 

>64 

>73 

>15 

Air  filled  with  foam  and  spray;  sea  completely  white 
with  driving  spray;  visibility  greatly  reduced 
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Table  7.  Saffir-Simpson  Scale  for  hurricane  intensity.  The  maximum  sustained  winds  are  given  in  meters  per  second,  knots,  and  miles 
per  hour. 


Category 

Sustained  winds 

Central 

pressure 

(mb) 

Storm 
surge  (m) 

Effects 

(ms  ^ 

(knots) 

(mph) 

Tropical 

storm 

17-32 

35-63 

39-73 

Beaufort  force  8-1 1 

1 

33-42 

64-82 

74-95 

>980 

1.0-1 .7 

No  real  damage  to  buildings.  Damage  primarily  to  unanchored  mobile  homes, 
shrubbery,  and  trees.  Some  flooding  of  coastal  roads  and  minor  damage  to  piers. 

2 

43 — 49 

83-95 

96-110 

979-965 

1. 8-2.6 

Some  damage  to  doors,  windows,  and  roofing  material.  Considerable  damage  to 
vegetation,  mobile  homes,  and  piers.  Coastal  low-lying  escape  routes  flood  2-4 
hours  before  the  storm  center  arrives.  Small  craft  in  unprotected  anchorages  break 
moorings. 

3 

50-58 

96-113 

111-130 

964-945 

2. 7-3.8 

Some  structural  damage  to  small  residences  and  utility  buildings,  with  a  minor 
amount  of  curtainwall  failures.  Mobile  homes  are  destroyed.  Flooding  near  the  coast 
destroys  smaller  structures,  with  larger  structures  damaged  by  floating  debris.  Terrain 
continuously  lower  than  5  feet  above  sea  level  may  be  flooded  8  miles  or  more 
inland. 

4 

59-69 

114-135 

131-155 

944-920 

3. 9-5. 6 

More  extensive  curtainwall  failures,  with  some  roofs  on  small  residences  failing 
completely.  Major  beach  erosion.  Major  damage  to  the  lower  floors  of  structures  near 
the  shore.  Terrain  continuously  lower  than  10  feet  above  sea  level  may  be  flooded 
inland  as  far  as  6  miles;  massive  evacuation  of  residential  areas  could,  therefore,  be 
required. 

5 

>69 

>135 

>155 

<920 

>5.7 

Complete  roof  failure  on  many  residences  and  industrial  buildings.  Some  complete 
building  failures,  with  small  utility  buildings  blown  over  or  away.  Major  damage  to 
lower  floors  of  all  structures  located  less  than  15  feet  above  sea  level  and  within  500 
yards  of  the  shore.  Massive  evacuation  may  be  required  for  residential  areas  on  low 
ground  within  5  to  1 0  miles  of  the  shore. 
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14.3  Fujita  Scale 

Similarly,  the  Fujita  Scale  categorizes  tornadoes  in  terms  of  their  maximum 
wind  speed  and  the  damage  they  cause  (Fujita  1981;  Glickman  2000,  p.  322f.). 
Table  8  shows  the  Fujita  Scale.  Fujita  (1981)  associated  the  lower  surface-level 
wind  speed  limit  for  a  tornado  category  with  the  Fujita  number  F  through  the 
expression 

U  =  6.30(F  +  2)3/2  ,  (14.2) 

where  U  is  the  wind  speed  in  m  s  '. 


Table  8.  Fujita  Scale  to  describe  tornado  intensity.  The  range  for  maximum  wind  speeds  is  given 
in  meters  per  second  and  miles  per  hour. 


Fujita 

scale 

Wind  speed  range 

Damage  specifications 

(m  s’1) 

(mph) 

F0 

18-32 

40-72 

Beaufort  force  8-1 1 .  Light  damage.  Some  damage  to  chimneys;  branches 
break  off  trees;  some  shallow-rooted  trees  pushed  over;  damage  to  sign 
boards. 

FI 

33 — 49 

73-112 

Moderate  damage.  The  lower  wind  speed  is  the  beginning  of  the  hurricane 
range.  Surfaces  of  roofs  peeled  off;  mobile  homes  pushed  off  foundations  or 
overturned;  moving  autos  pushed  off  the  road. 

F2 

50-69 

113-157 

Considerable  damage.  Roofs  torn  off  frame  houses;  mobile  homes 
demolished;  boxcars  pushed  over;  large  trees  snapped  off  or  uprooted;  light- 
object  missiles  generated. 

F3 

70-92 

158-206 

Severe  damage.  Roofs  and  some  walls  torn  off  well-constructed  houses;  trains 
overturned;  most  trees  in  a  forest  uprooted;  heavy  cars  lifted  off  the  ground 
and  thrown. 

F4 

93-116 

207-260 

Devastating  damage.  Well-constructed  houses  leveled;  structures  with  weak 
foundations  blown  off  some  distance;  cars  thrown;  large  missiles  generated. 

F5 

117-142 

261-318 

Incredible  damage.  Strong  frame  houses  lifted  off  foundations  and  carried 
considerable  distance  to  disintegrate;  automobile-sized  missiles  fly  through  the 
air  in  excess  of  100  m;  bark  ripped  from  trees;  incredible  phenomena  occur. 

F6+ 

>142 

>318 

Tornadoes  are  not  expected  to  reach  F6  wind  speeds. 
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15  CONCLUSIONS 

For  the  functions  that  I  have  presented  that  are  based  on  fits  to  data,  I  have 
generally  tried  to  estimate  the  uncertainty  in  the  predicted  values.  In  preparing 
this  review,  however,  I  was  surprised  to  see  how  old  some  of  the  data  are  that 
constitute  the  bases  for  these  fits.  For  example,  Flilsenrath  et  al.  (1960),  from 
which  1  take  several  functions,  was  published  more  than  40  years  ago.  Likewise, 
although  my  copy  of  the  Smithsonian  Meteorological  Tables  (List  1984)  has  a 
1984  publication  date,  this  volume  is  essentially  just  a  reprinting  of  the  original 
1951  tables.  Although  I  have  no  reason  to  doubt  the  quality  of  these  older  data, 
repeating  measurement  of  some  of  these  important  thermophysical  quantities 
with  modem  instruments  should  reduce  the  uncertainties  that  I  have  identified. 

In  particular,  many  of  the  gas  diffusivities  listed  in  Table  5  are  based  on  just 
theory  or  only  a  few  measurements.  The  gases  in  this  table — and  others  for  which 
I  could  find  no  information — are  environmentally  important  for  several  reasons: 
Some  are  greenhouse  gases,  some  are  pollutants,  and  others  are  useful  as 
environmental  tracers.  Lienee,  better  data  on  their  behaviors  in  air  is  crucial.  For 
example,  the  data  are  generally  so  scanty  that  we  still  cannot  decide  how  to 
predict  the  temperature  dependence  of  the  diffusivity,  as  reflected  in  the  variety 
of  exponents  in  (12.1),  (12.4),  and  (12.5).  In  my  view,  measuring  these 
diffusivities  over  a  range  of  temperatures  and  pressures  is  vital  research. 


Physical  Constants  and  Functions  for  Atmospheric  Boundary  Layer  Studies 


39 


REFERENCES 

Andreas,  E.L  (1989)  Thermal  and  size  evolution  of  sea  spray  droplets.  CRREL 
Report  89-11,  U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory, 
Hanover,  New  Hampshire,  37  p.  (NTIS:  AD  A2 10484.) 

Andreas,  E.L  (1990)  Time  constants  for  the  evolution  of  sea  spray  droplets. 

Tel l us,  42B:  481-497. 

Andreas,  E.L  (1995)  The  temperature  of  evaporating  sea  spray  droplets.  Journal 
of  the  Atmospheric  Sciences,  52:  852-862. 

Batchelor,  G.K.  (1970)  An  Introduction  to  Fluid  Dynamics.  Cambridge,  U.K.: 
Cambridge  University  Press,  615  p. 

Bohren,  C.F.,  and  B.A.  Albrecht  (1998)  Atmospheric  Thermodynamics.  New 
York:  Oxford  University  Press,  402  p. 

Bowditch,  N.  (1977)  American  Practical  Navigator,  Vol.  1.  Publication  No.  9, 
Washington,  D.C.:  Defense  Mapping  Agency  Hydrographic  Center,  1386  p. 

Brock,  F.V.,  and  S.J.  Richardson  (2001)  Meteorological  Measurement 
Systems.  New  York:  Oxford  University  Press,  290  p. 

Bromley,  L.A.,  V.A.  Desaussure,  J.C.  Clipp,  and  J.S.  Wright  (1967)  Heat 
capacities  of  sea  water  solutions  at  salinities  of  1  to  12%  and  temperatures  of  2° 
to  80°C.  Journal  of  Chemical  and  Engineering  Data,  12:  202-206. 

Buck,  A.L.  (1981)  New  equations  for  computing  vapor  pressure  and 
enhancement  factor.  Journal  of  Applied  Meteorology >,  20:  1527-1532. 

Businger,  J.A.  (1982)  The  fluxes  of  specific  enthalpy,  sensible  heat  and  latent 
heat  near  the  Earth’s  surface.  Journal  of  the  Atmospheric  Sciences,  39:  1 889— 
1892. 

Cox,  R.A.,  and  N.D.  Smith  (1959)  The  specific  heat  of  sea  water.  Proceedings 
of  the  Royal  Society,  London,  A252:  51-62. 

Fleagle,  R.G.,  and  J.A.  Businger  (1980)  An  Introduction  to  Atmospheric 
Physics.  2nd  ed.  New  York:  Academic  Press,  432  p. 

Fuehrer,  P.L.,  and  C.A.  Friehe  (2002)  Flux  corrections  revisited.  Boundary- 
Layer  Meteorology’,  102: 415-457. 

Fujita,  T.T.  (1981)  Tornadoes  and  downbursts  in  the  context  of  generalized 
planetary  scales.  Journal  of  the  Atmospheric  Sciences,  38:  1511-1534. 


40 


ERDC/CRREL  M-05-1 


Fuller,  E.N.,  K.  Ensley,  and  J.C.  Giddings  (1969)  Diffusion  of  halogenated 
hydrocarbons  in  helium.  The  effect  of  structure  on  collision  cross  section. 

Journal  of  Physical  Chemistry >,  73:  3679-3685. 

Gill,  A.E.  (1982)  Atmosphere-Ocean  Dynamics.  New  York:  Academic  Press, 

662  p. 

Glickman,  T.S.  (Ed.)  (2000)  Glossary  of  Meteorology’.  2nd  ed.  Boston: 

American  Meteorological  Society,  855  p. 

Goldstein,  S.  (Ed.)  (1965)  Modern  Developments  in  Fluid  Dynamics,  Vol.  1. 

New  York:  Dover,  330  p. 

Hall,  W.D.,  and  H.R.  Pruppacher  (1976)  The  survival  of  ice  particles  falling 
from  cirrus  clouds  in  subsaturated  air.  Journal  of  the  Atmospheric  Sciences,  33: 
1995-2006. 

Hanel,  G.  (1976)  The  properties  of  atmospheric  aerosol  particles  as  functions  of 
the  relative  humidity  at  thermodynamic  equilibrium  with  the  surrounding  moist 
air.  Advances  in  Geophysics,  19,  New  York:  Academic  Press,  73-188. 

Hare,  D.E.,  and  C.M.  Sorensen  (1987)  The  density  of  supercooled  water.  II. 
Bulk  samples  cooled  to  the  homogeneous  nucleation  limit.  Journal  of  Chemical 
Physics,  87:  4840-4845. 

Hilsenrath,  J,,  C.W.  Beckett,  W.S.  Benedict,  L.  Fano,  H.J.  Hoge,  J.F.  Masi, 
R.L.  Nuttall,  Y.S.  Touloukian,  and  H.W.  Woolley  (1960)  Tables  of 
Thermodynamic  Transport  Properties  of  Air,  Argon,  Carbon  Dioxide,  Carbon 
Monoxide,  Hydrogen,  Nitrogen,  Oxygen,  and  Steam.  Oxford,  U.K.:  Pergamon 
Press,  478  p. 

Hobbs,  P.V.  (1974)  Ice  Physics.  Oxford,  U.K.:  Clarendon  Press,  837  p. 

Horne,  R.A.  (1969)  Marine  Chemistry. \  New  York:  Wiley-lnterscience,  568  p. 

Iribarne,  J.V.,  and  W.L.  Godson  (1981)  Atmospheric  Thermodynamics.  2nd  ed. 
Dordrecht:  D.  Reidel,  259  p. 

Kester,  D.R.  (1974)  Comparison  of  recent  seawater  freezing  point  data.  Journal 
of  Geophysical  Research,  79:  4555-4556. 

Kinsman,  B.  (1965)  Wind  Waves.  Englewood  Cliffs,  New  Jersey:  Prentice-Hall, 
676  p. 

Larsen,  S.E.,  and  N.E.  Busch  (1974)  Hot-wire  measurements  in  the  atmosphere: 
Part  1.  Calibration  and  response  characteristics.  DISA  Information,  16:  15-34. 

Lide,  D.R.  (Ed.)  (2001)  CRC  Handbook  of  Chemistry  and  Physics.  82nd  ed. 

Boca  Raton,  Florida:  CRC  Press. 


Physical  Constants  and  Functions  for  Atmospheric  Boundary  Layer  Studies 


41 


List,  R.J.  (1984)  Smithsonian  Meteorological  Tables.  Sixth  ed.  Washington, 
D.C.:  Smithsonian  Institution  Press,  527  p. 

Lumley,  J,L,,  and  H.A.  Panofsky  (1964)  The  Structure  of  Atmospheric 
Turbulence.  New  York:  Wiley-lnterscience,  239  p. 

Massman,  W.J.  (1998)  A  review  of  the  molecular  diffusivities  of  H20,  C02, 
CH4,  CO,  O3,  S02,  NH3,  N20,  NO,  and  N02  in  air,  02,  and  N2  near  STP. 
Atmospheric  Environment,  32:  111  1-1 127. 

Millero,  F.J,,  G.  Perron,  and  J.E.  Desnoyes  (1973)  Heat  capacity  of  seawater 
solutions  from  5°  to  35°C  and  0.5  to  22%o  chlorinity.  Journal  of  Geophysical 
Research,  78:  4499-4507. 

Murphy,  D.M.,  and  T,  Koop  (2005)  Review  of  the  vapour  pressures  of  ice  and 
supercooled  water  for  atmospheric  applications.  Quarterly  Journal  of  the  Royal 
Meteorological  Society,  131:  1539-1565. 

Neumann,  G.,  and  W.J.  Pierson,  Jr.  (1966)  Principles  of  Physical 
Oceanography.  Englewood  Cliffs,  N.J.:  Prentice-Hall,  545  p. 

Pruppacher,  H.R.,  and  J.D.  Klett  (1978)  Microphysics  of  Clouds  and 
Precipitation.  Dordrecht:  D.  Reidel,  714  p. 

Pruppacher,  H.R.,  and  J.D.  Klett  (1997)  Microphysics  of  Clouds  and 
Precipitation.  2nd  revised  ed.  Dordrecht:  Kluwer,  954  p. 

Reid,  R.C.,  J.M.  Prausnitz,  and  B.E.  Poling  (1987)  The  Properties  of  Gases 
and  Liquids.  4th  ed.  New  York:  McGraw-Hill,  741  p. 

Reif,  F.  (1967)  Statistical  Physics:  Berkeley  Physics  Course,  Vol.  5.  New  York: 
McGraw-Hill,  398  p. 

Roll,  H.U.  (1965)  Physics  of  the  Marine  Atmosphere.  New  York:  Academic 
Press,  426  p. 

Schwerdtfeger,  P.  (1976)  Physical  Principles  of  Micro-Meteorological 
Measurements.  Amsterdam:  Elsevier,  113  p. 

Strangeways,  I.  (2001)  Back  to  basics:  The  ‘met  enclosure:’  Part  6-Wind. 
Weather,  56:  154-161. 

Thibodeaux,  L.J.  (1979)  Chemodynamics:  Environmental  Movement  of 
Chemicals  in  Air,  Water,  and  Soil.  New  York:  Wiley-lnterscience,  501  p. 

Thibodeaux,  L.J.  (1996)  Environmental  Chemodynamics:  Movement  of 
Chemicals  in  Air,  Water,  and  Soil.  2nd  ed.  New  York:  John  Wiley  and  Sons, 

593  p. 


42 


ERDC/CRREL  M-05-1 


Vargaftik,  N.B.,  B.N.  Volkov,  and  L.D.  Voljak  (1983)  International  tables  of 
the  surface  tension  of  water.  Journal  of  Physical  and  Chemical  Reference  Data, 
12: 817-820. 

Wagner,  P.E.  (1982)  Aerosol  growth  by  condensation.  Aerosol  Microphysics  II: 
Chemical  Physics  of  Microparticles  (W.H.  Marlow,  ed.).  Berlin:  Springer- 
Verlag,  129-178. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining 
the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for 
reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington, 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not 
display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YYYY) 
October  2005 

2.  REPORT  TYPE 

Monograph 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE 

5a.  CONTRACT  NUMBER 

Handbook  of  Physical  Constants  and  Functions  for  Use  in 
Atmospheric  Boundary  Layer  Studies 

_ 

_ 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

Edgar  L  Andreas 

_ 

_ 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Engineer  Research  and  Development  Center 

Cold  Regions  Research  and  Engineering  Laboratory 

72  Lyme  Road 

Hanover,  NH.  03755-1290 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

ERDC/CRREL  M-05-1 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Department  of  the  Army 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

Office  of  Naval  Research 

11. 

SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited 

Available  fromNTIS,  Springfield,  Virginia  22161 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

Studies  of  the  atmospheric  boundary  layer  always  require  values  for  the  dynamic  and  thermodynamic  properties  of  the  fluids  present: 
namely,  air,  water  vapor,  and  other  trace  gases.  Studies  of  the  marine  boundary  layer  frequently  require  similar  properties  for  the  ocean, 
in  addition.  This  handbook  collects  functions  for  calculating  the  values  for  these  dynamic  and  thermodynamic  variables.  For  example,  it 
includes  equations  for  air  density,  water  vapor  density,  and  seawater  density.  It  also  describes  various  humidity  variables  and  shows 
how  to  convert  among  these.  It  includes  equations  for  the  specific  heats  of  dry  air,  water  vapor,  moist  air,  seawater,  and  ice  and  for  the 
latent  heats  of  water  in  its  various  phases.  The  handbook  goes  on  to  show  how  to  calculate  the  molecular  transport  properties  such  as  the 
thermal  conductivity  of  air  and  the  diffusivities  of  heat,  water  vapor,  and  several  environmentally  important  gases  in  air.  This  discussion 
naturally  includes  values  for  the  molecular  Prandtl  and  Schmidt  numbers.  Finally,  the  handbook  concludes  with  summaries  of  several 
wind  force  scales  that  descriptively  categorize  winds  on  land  and  sea:  the  Beaufort  Scale,  the  Saffir-Simpson  Scale  for  hurricanes,  and 
the  Fujita  Scale  for  tornadoes. 

15.  SUBJECT  TERMS 

Air  density 

Beaufort  Scale 

Diffusivities  in  air  Mean  free  path  Surface  tension  of  water  Water  density 

Fujita  Scale  Saffir-Simpson  Scale  Thennal  conductivity  of  air  Water  vapor  variables 

Latent  heat  Specific  heat  Viscosity 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 
OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE 

PERSON 

a.  REPORT 

u 

b.  ABSTRACT 

u 

c.  THIS  PAGE 

u 

u 

52 

19b.  TELEPHONE  NUMBER  (include 
area  code) 

Standard  Form  298  (Rev.  8-98) 


Prescribed  by  ANSI  Std.  239.18 


